Variational density functional calculations of excited states via direct optimization

The development of variational density functional theory approaches to excited electronic states is impeded by limitations of the commonly used self-consistent field (SCF) procedure. A method based on a direct optimization approach as well as the maximum overlap method is presented and the performance compared with previously proposed SCF strategies. Excited-state solutions correspond to saddle points of the energy as a function of the electronic degrees of freedom. The approach presented here makes use of a preconditioner determined with the help of the maximum overlap method to guide the convergence on a target nth-order saddle point. The method is found to be more robust and to converge faster than previously proposed SCF approaches for a set of 89 excited states of molecules. A limited-memory formulation of the symmetric rank-one method for updating the inverse Hessian is found to give the best performance. A conical intersection for the carbon monoxide molecule is calculated without resorting to fractional occupation numbers. Calculations on excited states of the hydrogen atom and a doubly excited state of the dihydrogen molecule using a self-interaction corrected functional are presented. For these systems, the self-interaction correction is found to improve the accuracy of density functional calculations of excited states

Calculations of Excited Electronic States by Converging on Saddle Points Using Generalized Mode Following

Variational calculations of excited electronic states are carried out by finding saddle points on the surface that describes how the energy of the system varies as a function of the electronic degrees of freedom. This approach has several advantages over commonly used methods especially in the context of density functional calculations, as collapse to the ground state is avoided and yet, the orbitals are variationally optimized for the excited state. This optimization makes it possible to describe excitations with large charge transfer where calculations based on ground state orbitals are problematic, as in linear response time-dependent density functional theory. A generalized mode following method is presented where an $n^{\text{th}}$-order saddle point is found by inverting the components of the gradient in the direction of the eigenvectors of the $n$ lowest eigenvalues of the electronic Hessian matrix. This approach has the distinct advantage of following a chosen excited state through atomic configurations where the symmetry of the single determinant wave function is broken, as demonstrated in calculations of potential energy curves for nuclear motion in the ethylene and dihydrogen molecules. The method is implemented using a generalized Davidson algorithm and an exponential transformation for updating the orbitals within a generalized gradient approximation of the energy functional. Convergence is found to be more robust than for a direct optimization approach previously shown to outperform standard self-consistent field approaches, as illustrated here for charge transfer excitations in nitrobenzene and N-phenylpyrrole, involving calculations of $4^{\text{th}}$- and $6^{\text{th}}$-order saddle points, respectively. Finally, calculations of a diplatinum and silver complex are presented, illustrating the applicability of the method to excited state energy curves of large molecules.Comment: 57 pages, 12 figures, submitted to the Journal of Chemical Theory and Computatio

Direct Optimization Method for Variational Excited-State Density Functional Calculations Using Real Space Grid or Plane Waves

A direct optimization method is presented for density functional calculations of excited electronic states using either a real space grid or a plane wave basis set. The method is variational, provides atomic forces in the excited states, and can be applied to Kohn-Sham (KS) functionals as well as orbital-density dependent functionals (ODD) including explicit self-interaction correction. The implementation for KS functionals involves two nested loops: (1) An inner loop for finding a stationary point in a subspace spanned by the occupied and a few virtual orbitals corresponding to the excited state; (2) an outer loop for minimizing the energy in a tangential direction. For ODD functionals, a third loop is used to find the unitary transformation that minimizes the energy functional among occupied orbitals only. Combined with the maximum overlap method, the algorithm converges in challenging cases where conventional self-consistent field algorithms tend to fail. The benchmark tests presented include two charge-transfer excitations in nitrobenzene and an excitation of CO to degenerate $\pi^\ast$ orbitals where the importance of complex orbitals is illustrated. An application of the method to several metal-to-ligand charge-transfer and metal-centred excited states of an Fe$^{\rm II}$ photosensitizer complex is described and the results compared to reported experimental estimates. The method is also used to study the effect of Perdew-Zunger self-interaction correction on valence and Rydberg excited states of several molecules, both singlet and triplet states. The correction is found to improve the description of molecular bond stretching but calculated values of the excitation energy are improved only slightly, by {\it ca.} 0.1 eV, due to cancellation of the estimated self-interaction error in the ground and excited states.Comment: 55 pages, 12 figures, including supporting Informatio

On how mild oxidation affects the structure of carbons: Comparative analysis by different techniques

Abstract Understanding how the structure of carbonaceous materials changes upon oxidation at mild temperature as a function of the graphitic order is of great importance for the development of clean combustion technologies, such as carbon fuel cells. The micro- and nanostructures of a range of carbon materials at room temperature and upon mild oxidation at 300 °C have been analysed by means of mercury porosimetry, Nitrogen adsorption, X-Ray Diffraction, Raman spectroscopy and Electron Paramagnetic Resonance. The samples included carbons with increasing level of graphitic order: three chars derived from two bituminous coals and a lignite, a synthetic carbon and a graphitized coke. The experimental characterization allows to classify the materials according to different structural parameters, including porosity, surface area, degree of graphitization and paramagnetic activation of the carbonaceous surface. A correlation with the rank of the analysed materials is observed. For the graphitized coke, oxidation leads to more crystalline order and enhancement of the paramagnetic signal. A similar increase in the paramagnetic activity is observed for the Lignite char. On the other hand, for the higher rank, bituminous and synthetic chars, mild oxidation leads to a slight expansion of the amorphous carbon and loss of paramagnetic activity. The differences are rationalized in terms of formation of new carbon-oxygen complexes on the graphitized coke and on the Lignite char, and redistribution of already existing complexes for the higher-rank coals. This investigation complements previous X-ray photoelectron spectroscopy measurements

Variational calculations of excited states via direct optimization of the orbitals in DFT

Post-print (lokagerð höfundar)A direct optimization method for obtaining excited electronic states using density functionals is presented. It involves selective convergence on saddle points on the energy surface representing the variation of the energy as a function of the electronic degrees of freedom, thereby avoiding convergence to a minimum and corresponding variational collapse to the ground electronic state. The method is based on an exponential transformation of the molecular orbitals, making it possible to use efficient quasi-Newton optimization approaches. Direct convergence on a target nth-order saddle point is guided by an appropriate preconditioner for the optimization as well as the maximum overlap method. Results of benchmark calculations of 52 excited states of molecules indicate that the method is more robust than a standard self-consistent field (SCF) approach especially when degenerate or quasi-degenerate orbitals are involved. The method can overcome challenges arising from rearrangement of closely spaced orbitals in a charge-transfer excitation of the nitrobenzene molecule, a case where the SCF fails to converge. The formulation of the method is general and can be applied to non-unitary invariant functionals, such as self-interaction corrected functionals.The present work was funded by the Icelandic Research Fund (grant number 196070-052) and the University of Iceland Research Fund. AVI is supported by a doctoral fellowship from the University of Iceland.Peer Reviewe

Excited-State Solvation Structure of Transition Metal Complexes from Molecular Dynamics Simulations and Assessment of Partial Atomic Charge Methods

In this work, we investigate the excited-state solute and solvation structure of $\mathrm{[Ru(bpy)_3]^{2+}}$, $\mathrm{[Fe(bpy)_3]^{2+}}$, $\mathrm{[Fe(bmip)_2]^{2+}}$ and $\mathrm{[Cu(phen)_2]^{+}}$ (bpy=2,2'-pyridine; bmip=2,6-bis(3-methyl-imidazole-1-ylidine)-pyridine; phen=1,10-phenanthroline) transition metal complexes (TMCs) in terms of solute-solvent radial distribution functions (RDFs) and evaluate the performance of some of the most popular partial atomic charge (PAC) methods for obtaining these RDFs by molecular dynamics (MD) simulations. To this end, we compare classical MD of a frozen solute in water and acetonitrile (ACN) with quantum mechanics/molecular mechanics Born-Oppenheimer molecular dynamics (QM/MM BOMD) simulations. The calculated RDFs show that the choice of a suitable PAC method is dependent on the coordination number of the metal, denticity of the ligands, and type of solvent. It is found that this selection is less sensitive for water than ACN. Furthermore, a careful choice of the PAC method should be considered for TMCs that exhibit a free direct coordination site, such as $\mathrm{[Cu(phen)_2]^{+}}$. The results of this work show that fast classical MD simulations with ChelpG/RESP or CM5 PACs can produce RDFs close to those obtained by QM/MM MD and thus, provide reliable solvation structures of TMCs to be used, e.g. in the analysis of scattering data

Orbital-optimized Density Functional Calculations of Molecular Rydberg Excited States with Real Space Grid Representation and Self-Interaction Correction

Density functional calculations of Rydberg excited states up to high energy are carried out for several molecules using an approach where the orbitals are variationally optimized by converging on saddle points on the electronic energy surface within a real space grid representation. Remarkably good agreement with experimental estimates of the excitation energy is obtained using the generalized gradient approximation (GGA) functional of Perdew, Burke and Ernzerhof (PBE) when Perdew-Zunger self-interaction correction is applied in combination with complex-valued orbitals. Even without the correction, the PBE functional gives quite good results despite the fact that corresponding Rydberg virtual orbitals have positive energy in the ground state calculation. Results obtained using the TPSS and r2SCAN meta-GGA functionals are also presented, but they do not provide a systematic improvement over the results from the uncorrected PBE functional. The grid representation combined with the projector augmented-wave approach gives a simpler and better representation of the diffuse Rydberg orbitals than a linear combination of atomic orbitals with commonly used basis sets, the latter leading to an overestimation of the excitation energy due to confinement of the excited states.Comment: 34 pages, 8 figure

Ethical and Social Challenges with developing Automated Methods to Detect and Warn potential victims of Mass-marketing Fraud (MMF)

Mass-marketing frauds (MMFs) are on the increase. Given the amount of monies lost and the psychological impact of MMFs there is an urgent need to develop new and effective methods to prevent more of these crimes. This paper reports the early planning of automated methods our interdisciplinary team are developing to prevent and detect MMF. Importantly, the paper presents the ethical and social constraints involved in such a model and suggests concerns others might also consider when developing automated systems